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Abstract. - Ideal bosons and classical ring polymers formed via self-assembly, are known 
to have the same partition function, and so analogous phase transitions. In ring polymers, 
the analogue of Bose-Einstein condensation occurs when a ring polymer of macroscopic size 
appears. We show that a transition of the same general form occurs within a whole class of 
systems with self-assembly, and illustrate it with the emulsification failure of a microemulsion 
phase of water, oil and surfactant. As with Bose-Einstein condensation, the transition occurs 
even in the absence of interactions. 



Bose-Einstein condensation (BEC) is a textbook pi, but rather unusual, phase transition. 
It occurs for noninteracting bosons, in contrast to more conventional transitions such as that 
of the Ising model, which are driven by interactions. Feynman 0] showed that the statistical 
mechanics of bosons can be performed via what he called path integrals. These integrals 
are, in turn, equivalent to integrals over the configurations of ring polymers. Thus, the 
path integral formalism of Feynman implies that noninteracting bosons and self-assembling, 
noninteracting ring polymers have partition functions of exactly the same form. Necessarily 
then, ring polymers must undergo a phase transition precisely analogous to BEC p]. Here, we 
generalise this result to show that there is a class of self-assembling systems which undergo a 
phase transition analogous to BEC. This phase transition occurs in the absence of interactions 
between the aggregates formed by self-assembly. At BEC a condensate appears which is a 
macroscopic number of bosons in a single state. In the analogous transition in self-assembling 
systems, an aggregate of macroscopic size appears. For the example considered here, a 
microemulsion, the macroscopic aggregate is a bulk oil phase. 

Self-assembling systems are systems in which the particles are not immutable objects but 
are formed reversibly Q Typically we have water, surfactant and sometimes oil. The 
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surfactant molecules then spontaneously assemble into micelles, or coat and stabilise droplets 
of oil. Both the micelles and the droplets are what we term aggregates. For example, in a 
microemulsion, equilibrium is obtained when the oil is dispersed in water as oil droplets whose 
surfaces are coated with surfactant. These droplets are more stable than just a single bulk 
oil phase because with a single bulk phase there is no extensive oil-water interface, and the 
amphiphilic surfactant molecules have the lowest energy at this interface. The distribution 
of droplets changes with density and temperature. A BEC-like transition occurs when, at 
equilibrium, some of the oil exists as a macroscopic droplet, i.e., a bulk phase. This transition 
is called emulsiftcation failure . 

Let us begin by outlining BEC in Feynman's language The partition function of N 
bosons is given in terms of their density matrix p by 

Z N (p,V)=tv v p, p = e- pHN , (1) 

with Hn the Hamiltonian operator, /3 = 1/kT. The density matrix of noninteracting identical 
bosons can be written in terms of 1-particle density matrices, /3i(xj;xQ, 

1 N 

(3 ({x i };{^}) = M J2 Il^( x - Px P(*))< ( 2 ) 
' Pen N i=i 

where {x^} = {xi, . . . , xjv} and II at is the set of permutations of iV elements [P : i i— * P(i)}. 

Taking the trace over eq. (||) and moving over to the grand canonical ensemble, we obtain 
for the grand partition function, S(/3, z, V), and average density of bosons, (j), 



ms = ( 3 ) 

n— 1 n— 1 

with z the activity of the bosons and h n = try = Zi(n(3,V). It is straightforward to 
determine that for a <i-dimensional cubic box with periodic boundary condition Zi(/3,V) = 
[$(e _7r ^ )] d , where £ = V x l d / A is the length of the box in units of the thermal wavelength 
A = h v / (3/2nm, and ti(q) = 1 + 2 Y^=i Q k * • Thus > 



1 



z n ^( e --™/^ 2 )] . ( 4 ) 



Remarkably, eq. (^) also applies to self-assembling ring polymers, provided z and A are 
properly redefined. To see this, note that by definition H, the one-particle density matrix, 
pi = exp(—[3Hi), is a solution of the equation dppi = —H\p\ = (ti 2 / '2m)V 2 p\ . This is just a 
Fokker-Planck equation for a random walk in "time" (3, and with diffusion coefficient h 2 /2m. 
Accordingly, l/ _1 try P\{— h n /nV) has an alternative interpretation as the probability of a 
random walk forming a cyclic path after n "time" -steps of length (3. The simplest model of a 
ring polymer of length n is a cyclic random walk of n steps; hence h n /nV can also be regarded 
as the ratio of the partition function of a ring of n monomers, Z" ns , to that of a chain of n + 1 
monomers, VA nd . With this interpretation A has the meaning of the average monomer length. 
We conclude that Z™ s = A nd h n /n. For an ideal mixture of rings made up of monomers at 
an activity z, the number density of rings of length n, p{n) — z n Z^ ng /V, and the number 
density of monomers in rings of length n is just n times this, (zA d ) n h n /V. If we absorb a 
factor of A d into z then clearly the total number density of monomers in ring polymers is 
given exactly by the expression for the density of bosons in eq. ([|) . This well-known mapping 
of ring polymers onto bosons is called the "classical isomorphism" ||. It tells us that both 
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systems have completely equivalent equilibrium behaviour, and so (i) the ring polymer system 
must undergo the analogue of BEC, and (ii) ideal bosons are equivalent to a classical system 
with self-assembling aggregates. 

To locate the phase transition for bosons/ring polymers let us consider an integer n\(V) 
such that m < I 2 but n x {V) -^ooasy^oo. Then i?(e" 7r "/^) = In- 1 ! 2 + 0(e-* e ) for all 
n < ni(V), and so we can split the sum of eq. (^), yielding 

1 ni n 1 00 d 

^A^Ei4^ Zn {^ n/e2 )} + <V), (5) 

n— 1 n—ni 

where e(V) — > as V — > oo. If we fix z < 1 and take the thermodynamic limit then the 
first sum is less than its value at z = 1, <f) c = A _d ^^° =1 n _d / 2 , and the second sum is zero. 
When z > 1 the first sum is clearly divergent. This implies that all cj> > <fi c correspond to 
z = 1. Careful analysis of eq. (||) shows || that for > <f> c , the density in excess of (f> c , 
i.e., 4> — <f> c , comes from the second sum of eq. (|^). This second sum is the contribution of 
either a macroscopic number of bosons in a single state (the ground state) or a ring polymer 
of macroscopic size. 

We now generalise the preceding theory to describe a whole class of systems with self- 
assembly. Let us consider monomers which can self-assemble into aggregates, where these 
aggregates are noninteracting but are otherwise arbitrary. The internal partition function of 
the corresponding aggregates in the new model will differ from that of ring polymers, but 
density of matter, </>, will again be given in terms of the activity, z, through some function 
4> = G(z). We will then have a BEC-like transition provided that as z runs over a range 
< z < z c , the new function takes only values in the limited range < G(z) < G(z c ), and is 
divergent for z > z c . Furthermore, this phase transition will again occur due to the appearance 
of macroscopic aggregates. 

Let us denote by s the (dimensionless) "size" (volume, surface, length ... ) of aggregates, 
and assume sq < s < oo. If p(s) denotes the number density of aggregates of size s, then the 
free energy of an ideal mixture of these aggregates is 

/3F/V = ^2p(s)[lnp(s)-l + f(s) + as] 

s>s 

+Popnpo-l + /o] (6) 

where f(s) is the internal free energy of aggregates of size s, with the linear part, as, subtracted 
off. The second term must be included if aggregates for which s — are present (for example, 
in microemulsions micelles are present which are pure surfactant and do not contribute to 
<j>; see below). In this term po and fo are the number density and internal free energy of 
these aggregates. The sum in eq. (JsJ) may be replaced by an integral with only very minor 
quantitative changes to the behaviour. We will do so when we come to do explicit calculations 
for microemulsions, for which s may be treated as a continuous variable. Note that ring 
polymers conform to eq. (|^) if s denotes the number of monomers (i.e. the dimensionless 
length), So = 1 and /(s) + as = — \n[Z rins (s)/V] (the s = term is obviously absent). 

A microemulsion is an isotropic equilibrium phase composed of water, oil and a surfactant 
@ §. If there is, say, much less oil than water then the microemulsion consists of droplets 
of oil, nanometers across, coated with the surfactant and dispersed in the water. Without 
the surfactant, at equilibrium we would have bulk oil and water phases, but as surfactant 
prefers to lie in the water-oil interface, it stabilises the oil droplets dispersed in the water. 
For microemulsions the sum of eq. (^) is over droplets of different sizes; here s is the reduced 
volume of a droplet and So is the size of the smallest possible droplet. Micelles, aggregates 
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of pure surfactant, are also present in microcmulsions, and as they are pure surfactant they 
do not contribute to the oil density <f>. In eq. ([3j) applied to microemulsions, the second term 
on the right hand side accounts for micelles; po is the density of micelles. The surfactant is 
assumed to be restricted to the surface of the droplets and to micelles. 

At equilibrium the droplet distribution p(s) is given by the minimum of the free energy, 
subject to the constraints that the total fraction of oil, 

</> = E s ^ s )' ( ? ) 

s>so 

is fixed, and also that the total amount of surfactant is fixed. Fixing the (reduced) surfactant 
density £ requires constraining the sum of the micelle density, po, and the s 2 ' 3 moment of p(s). 
The latter is proportional to the surface area of the droplets per unit volume. This constraint 
is a single specific example of constrained densities £j of the form 

& = c^po + ^2 w l (s)p(s), i = l,...,v, (8) 

s>s 

where Cj are constants, po is a density of a species which contributes to & but not to and 
Wi(s) are weight functions. Equations (^), ([?]) and (|J) define our model of a system with 
self-assembly. 

Minimising the general free energy, eq. (||), at constant </>, eq. (|?|), under the constraints (||), 
we obtain the equilibrium densities 

p(s) = e -/w -£« A ^*( s )- As , pa = e -A-E« x ^ , (9) 

where X — a — ftp, p being the chemical potential corresponding to the reduced density <j>, and 
Xi are the Lagrange multipliers corresponding to the constraints (H) (and which are analogous 
to chemical potentials associated with the fixed densities 

Now, if either /(s) or any of the u>j(s) grows faster than s as s — > oo, then the sum in eq. (|7]) 
will converge for any value of A. Furthermore, for any set of densities {</>, £j} there will always 
be a set {A, Ai}, which will solve eqs. (^) and (]?]). This means that the free energy will be a 
smooth function of the densities and no phase transition will occur. Therefore, for a phase 
transition to occur, it is necessary that (a) f(s), Wi(s) — o(s) as s — > oo. If this condition holds 
then the sum in eq. (|?j) is convergent for A > and divergent for A < 0. The second condition 
concerns the limit A — > + . If the sum in eq. ([?]) diverges in this limit, then again any set 
of densities will be reachable and no phase transition will occur. Thus, for the transition to 
appear it is necessary that (b) the sum in eq. (0) is convergent for A = 0. 

If our model fulfills conditions (a) and (b) we have that the density <f> increases to a finite 
value and then suddenly diverges, just as for bosons. Just as with bosons, to study the 
transition we need to consider explicitly the volume dependence. For this we need to consider 
a free energy of aggregates which depends on V, i.e., lnp(s) — l+/(s; V)+as. Now, f(s;V) has 
to be such that /(s; V) — > f(s) as V — ► oo, but the convergence cannot be uniform: intuitively, 
if s is "small" compared to V, then f(s; V) « /(s), but if s is "large", then /(s; V) will differ 
markedly from f(s). Thus we can always choose a size si(V) such that f(s < si;V) ~ f(s) 
and that s±(V) — > oo when V — * oo. 

There are two possible scenarios depending on whether /(s; V) violates either condition (a) 
or condition (b). For example, condition (a) is violated when compact aggregates form and 
the box limits the maximum size of an aggregate. This translates into a sharp increase of 
/(s; V) at a certain macroscopic size. Case (b) is what we have met for ring polymers, for 
which f(s; V) ~ ln(Vs) as s — > oo. Although there are small differences between the two cases 
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[Bj, they do not affect the dominant terms in the thermodynamic limit. We report the details 
of the analysis elsewhere ||, and limit ourselves to describing the resulting phase behaviour. 
If <j) < cf) c , eqs. (0) and (||) can be solved for A and A^. This yields the densities, eq. (^|), as well 
as the free energy 

/3F/V=-p-J2 X ^i+ a ^ ( 10 ) 

i=l 

with p = po + y^ s p( s )- The phase transition occurs at a density <p c , which is given by setting 
A = in cq. (g) for p and inserting the result in eq. (0). When cf> > cf> c , the matter in excess 
of 4>ci 4> ~ 4>c forms macroscopically big aggregates (the condensate), and A = 0. The density 
of non- macroscopic aggregates is then given by eq. (^|) with A = and A.; = Ai. c , where 
Xi.c are the solution to eqs. (||) for A = 0. In particular, as /(s) increases sublinearly then 
whenever macroscopic aggregates are present, A = and the size distribution of aggregates 
always decays slower than exponentially. Finally, in this regime the free energy has the same 
form as in eq. (|i"o|), but with p and Ai replaced by p c and Xi. c , the corresponding values for 
A = 0. 

Returning to our example of microemulsions, we will explicitly calculate a phase diagram. 
As mentioned before, s is proportional to the volume. The internal free energy of a droplet has 
a linear term, as, plus /(s), which comes from the droplet's surface. The surface of the drop 
can be modeled using a surface tension term plus an elastic free energy of the form introduced 
by Helfrich Jll], which accounts for a preferred curvature of the surfactant layer at the 
surface of the droplet. The details are given in Ref. ||, but as we might expect there are 
three contributions arising from the surface tension (proportional to s 2 / 3 ), the mean curvature 
(proportional to s 1 / 3 ) and the Gaussian curvature (a constant ^ 1 )[ Thus we take for /(s) 

f{s)^a -a 1 s 1 / 3 +a 2 s 2 / 3 , (11) 

consistent with requirement (a) above; a\ > 0, the signs of ao and a 2 are not fixed. This 
expression for f(s) may be inserted in the expression for the free energy, eq. (|^), and the sum 
replaced by an integral; s may be any volume in the range sq < s < 00. There is a single 
constraint, that on the surfactant density 

roc 

e = c po+ci / dss 2/3 p(s) (12) 

J s 

(co and c\ are simply geometric factors). Minimising the free energy subject to fixed <f> and 
£ yields the droplet size distribution p(s) = cxp{— a + ais 1 / 3 — (a 2 + A 2 ci)s 2 / 3 — As}, and 
density of micelles po — e~^°~ c ° A2 . If a 2 + A 2 ci > the requirement (b) on /(s) also holds 
and we will have a finite <p c at which a transition occurs. At (j) c a macroscopic droplet forms, 
which is nothing other than a bulk oil phase. So, a bulk phase of excess oil coexisting with 
the microemulsion phase has formed. This is the so-called emulsification failure 0. 

To illustrate emulsification failure we take some simple and rather arbitrary values for the 
parameters of the free energy, namely a = ao = 0,2 = fo = 0, a\ = c$ = c\ = sq = 1, 
and calculate the phase diagram, fig. [l]. At fixed surfactant density £ there is a maximum 
oil density C beyond which the microemulsion phase coexists with a bulk oil phase. This 
maximum density increases as the surfactant density increases. 

( 1 ) Reiss and coworkers jl2| have looked carefully at the function f(s) for microemulsions and 
found in addition to the terms in eq. ( pd[ ) a term logarithmic in s. This term essentially accounts for 
fluctuations of the center of mass of thearoplet. Adding such a term to eq. (fill) results in a /(s) which 
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Fig. 1. - Amount of oil in droplets, <f>, vs. amount of surfactant, f. In the shaded region a bulk oil 
phase coexists with the microemulsion (emulsification failure). The solid curve (A = 0) signals the 
phase transition. The dashed curves are at constant micelle density po = e~ A2 . From bottom to top 
A 2 = -1,0,1,2. 



Finally we briefly comment on intcraggregate interactions. Their presence will of course 
affect any transition which is driven by the mechanism outlined in this Letter. But in 
addition, Zhang et al. fl3|| report evidence of a transition like that found here for a system 
for which f(s) =const.. There, the interactions between the aggregates lead to a distribution 
of aggregates of the form of eq. (||) but with f(s) depending on the density as a result of the 
interactions, not just arising from an internal free energy [[l3| 14, [l^]. If this f(s) satisfies 



conditions (a) and (b) there will be a transition induced by the interactions. 

There is a class of systems in which aggregates of all sizes can reversibly self-assemble, 
and in which, without interaggregate interactions, there is a phase transition. This phase 
transition occurs when an infinite aggregate or aggregates form. It is analogous to the 
Bose-Einstein condensation of ideal bosons: in neither case are there interactions, and the 
macroscopic aggregates we find play the same role as the condensate of a macroscopic number 
of bosons. This analogy for the specific case when the aggregates are ring polymers is implied by 
Feynman's path integral approach for bosons [|[ ||, H) . We have shown that the analogy applies 
to a whole class of self-assembling systems with partition functions which are qualitatively 
similar, although of a different functional form, to that of ideal bosons/ring polymers. This 
class is defined by the requirement that the internal free energy of an aggregate contains a 
linear term, as, plus others which are sublinear and ensure that the density is finite when 
the chemical potential /3/i — > a~ . We outlined the behaviour of one example here, that of 



still satisfies condition (a) for a transition, and thus does not alter the nature of the phase behaviour. 
It will however shift the phase boundary. 
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microemulsions. Another example is that of the formation of the lamellar phase, a stack 
of effectively infinite bilayers of surfactant, which can coexist with small discs of bilayer Jl6[ . 
This transition, as well as emulsification failure, will be affected by the presence of interactions 
between the aggregates. We suggest that this effect will be analogous to that of interactions 
on the A transition of liquid 4 Hc pL thought to be Bose- Einstein condensation modified by 
the presence of interactions between the He atoms. 
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